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Abstract. We study the behavior of steady state voltage potentials in two 
kinds of bidimensional media composed of material of complex permittivity 
equal to 1 (respectively a) surrounded by a thin membrane of thickness h and 
of complex permittivity a (respectively 1). We provide in both cases a rigor- 
ous derivation of the asymptotic expansion of steady state voltage potentials 
at any order as h tends to zero, when Neumann boundary condition is im- 
posed on the exterior boundary of the thin layer. Our complex parameter a 
is bounded but may be very small compared to 1, hence our results describe 
the asymptotics of steady state voltage potentials in all heterogeneous and 
highly heterogeneous media with thin layer. The terms of the potential in the 
membrane are given explicitly in local coordinates in terms of the boundary 
data and of the curvature of the domain, while these of the inner potential are 
the solutions to the so-called dielectric formulation with appropriate bound- 
ary conditions. The error estimates are given explicitly in terms of h and a 
with appropriate Sobolev norm of the boundary data. We show that the two 
situations described above lead to completely different asymptotic behaviors 
of the potentials. 



Introduction 

We study the behavior of the steady state voltage potentials in highly contrasted 
media surrounded by a thin layer. The motivation of the present work comes from 
numerical problems raised by the researchers in computational electromagnetics, 
who want to compute the quasi-static electric field in highly contrasted materials 
with thin layer. The thinness of the membrane surrounding an inner domain leads 
to numerical difficulties, in particular for the meshing. 

To avoid these difficulties, we perform an asymptotic expansion of the potentials 
in terms of the membrane thickness. The approached inner potential is then the 
finite sum of the solutions to elementary problems in the inner domain with appro- 
priate conditions on its boundary, which approximate the effect of the thin layer. 
Thereby, the thin membrane does not have to be considered anymore. Our method 
leads to the construction of so-called "approximated boundary conditions" at any 
order [11] . We estimate precisely the error performed by this method in terms of 
an appropriate power of the relative thinness and with a precise Sobolev norm of 
the boundary data. This method is well-known for non highly contrasted media. 
It is formally described in some particular cases in p] and [15] . We also refer to 
Krahenbuhl and Muller [TH] for electromagnetic considerations. Usually, when it is 
estimated (see for example |11| ) , the norm of the error involves a imprecise norm 
of the boundary data (a < ^' 00 norm while a weaker norm is enough) and mainly, 
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the constant of the estimate depends strongly on the dielectric parameters of the 
domain. 

It is not obvious (it is even false in general!) that such results hold for highly 
contrasted domains with thin layer, and this is a fact that researchers in computa- 
tional electromagnetics are often confronted to such media. For example, a simple 
electric modelization of the biological cell consists of a conducting cytoplasm sur- 
rounded by a thin insulating membrancQ; the modulus of the cytoplasmic complex 
permittivity divided by the membrane permittivity is around 10 5 while the relative 
thinness is equal to 10~ 3 . On the other hand, the medium might be a dielectric 
surrounded by a thin metallic layer. In both cases it is not clear that the usual 
approximated boundary condition might be used. 

We derive asymptotics of the potential steady state voltage in all possible do- 
mains with thin layer (heterogeneous or highly heterogeneous). As we will see, the 
two situations described above lead to different behaviors of the potentials. The 
membrane relative thickness is equal to h, while the charateristic length of the in- 
ner domains is equal to 1. The first medium consists of a conducting inner domain 
(say that its complex permittivity is equal to 1) surrounded by a thin membrane; 
we denote by a the membrane complex permittivity. The parameter a is bounded 
but it may tend to zero. This is the reason why we say that the thin layer is an 
insulating membrane. The second material consists of an insulating inner domain 
of permittivity a surrounded by a conducting thin membrane (say that its complex 
permittivity is equal to 1). In this case, we suppose that a tends to zero. These 
two kinds of media describe all the possible media with thin layer. The aim of 
this paper is to derive full rigorous asymptotic expansion of steady state voltage 
potentials with respect to the small parameter h for bounded a (but it may tend 
to zero). 

Let us write mathematically our problem. Let flh be a smooth bounded bidi- 
mensional domain (see Fig. |T|) , composed of a smooth domain O surrounded by a 
thin membrane Oh with a small constant thickness h: 



Let a be a non null complex parameter with positive real part; a is bounded but 
it may be very small. Without loss of generality, we suppose that \a\ < 1. Denote 
by qh and "/h the following piecewise constant functions 



We would like to understand the behavior for h tending to zero and uniformly 
with respect to |a| < 1 of Vh and ut the respective solutions to the following 
problems ([T]) and ((2J) with Neumann boundary condition; Vh satisfies 



n h = ouo h . 




(la) 



V • (q h VV h ) 

dV h 
dn 



in Q h , 



(lb) 



(j) on dilh, 



(lc) 




0: 



We refer to the author thesis 1171 for a precise description of the biological cell. 
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and Uh satisfies 

(2a) V • { lh Vu h ) = in Q h , 
(2b) ^ = on dCl h , 

(2c) / u h da = 0. 

JdO 

Since we impose a Neumann boundary conditions on dQh the boundary data <fi 
must satisfy the compatibility condition: 

/ 0dc = O. 

The above functions Vh and Uh are well-defined and belong to if 1 (ft/,) as soon as 
<j> belongs to H- 1 / 2 {dVL h ). 

Several authors have worked on similar problems (see for instance Beretta et 
al. [5] and [6]). They compared the exact solution to the so-called background 
solution defined by replacing the material of the membrane by the inner material. 
The difference between these two solutions has then been given through an integral 
involving the polarization tensor defined for instance in [5] , [3] , [5] , [S] , ■> phis some 
remainder terms. The remainder terms are estimated in terms of the measure of 
the inhomogeneity. In this paper, we do not use this approach, for several reasons. 

The Beretta et al. estimate of the remainder terms depends linearly on a and 
I /a: their results are no more valid in a highly contrasted domain (i.e. for a very 
large or very small). Secondly, a is complex- valued, hence differential operators 
involved in our case are not self-adjoint, so that the T-convergence techniques of 
Beretta et al. do not apply. Thirdly, the potential in the membrane is not given 
explicitly in [5], [6] or [7], while we are definitely interested in this potential, in 
order to obtain the transmembranar potential (see Fear and Stuchly [H]). Finally, 
the asymptotics of Beretta et al. are valid on the boundary of the domain, while 
we are interested in the potentials in the inner domain. 

The heuristics of this work consist in performing a change of coordinates in the 
membrane Oh, so as to parameterize it by local coordinates (f),0), which vary in 
a domain independently of h; in particular, if we denote by L the length of dO 
(in the following, without any restriction, we suppose that L is equal to 2-7r)), the 
variables (77, 9) should vary in [0, 1] x M./LZ. This change of coordinates leads to an 
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expression of the Laplacian in the membrane, which depends on h. Once the trans- 
mission conditions of the new problem are derived, we perform a formal asymptotic 
expansion of the solution to Problem (p} (respectively to Problem (J3J)) in terms of 
h. It remains to validate this expansion. In this paper we work with bidimcnsional 
domain and we are confident that the same analysis could be perfomed in higher 
dimensions. 

This paper is structured as follows. In Section [TJ we make precise our geometric 
conventions. We perform a change of variables in the membrane, and with the help 
of some differential geometry results, we write Problem (JT]) and Problem ([2|) in the 
language of differential forms. We refer the reader to Flanders [13] or Dubrovin 
et al. [5] (or [5] for the french version) for courses on differential geometry. We 
derive transmission and boundary conditions in the intrinsic language of differential 
forms, and we express these relations in local coordinates. 

In Section [3] we study Problem In paragraph 13. II we derive formally all the 
terms of the asymptotic expansion of the solution to our problem in terms of h. 
Paragraph 13. 21 is devoted to a proof of the estimate of the error. 

Problem ^ is considered in Section HJ We supposed that a tends to zero: a 
boundary layer phenomenon appears. To obtain our error estimates, we link the 
parameters h and a. We introduce a complex parameter [3 such that 

Re(/3) > 0, or (Re(/3) = 0, and 3(/3) ^ 0) , 

and 

We distinguish two different cases, depending on the convergence of |a| to zero: 
a = [3h q , for q 6 N* and a — o(h N ) for all N e N. 

For g=lwe obtain mixed boundary conditions for the asymptotic terms of the 
inner potential, and as soon as q > 1, appropriate Dirichlet boundary conditions 
are obtained. We end this section by error estimates. 

In Appendix, we give some useful differential geometry formulae. 

Remark 1. The use of the formalism of differential forms 5 ((fad) could seem futile 
for the study of the operator V • ((faV). In particular the expression of Laplace 
operator in local coordinates is well known. However we wanted to present this 
point of view to show how simple it is to write a Laplacian in curved coordinates 
once the metric is known. 

Moreover once this formalism is understood for the functions ( or 0-forms ), it 
is easy to study 5 ((fad) applied to 1-forms. This leads directly to the study of the 
operator rot ((fa rot) , whose expression in local coordinates is less usual. 

We choose to present our two main theorems in this introduction so that the 
reader interested in our results without their proves might find them easily. 

We suppose that dO is smooth. We denote by <E> the C € (X '— diffeormorphism, 
which maps a neighborhood of cylinder C = [0, 1] x M/27rZ unto a neighborhood of 
the thin layer. The diffeomorphism (f>o — ^(0, •) maps the torus unto the boundary 
dO of the inner domain while $i = $(1-) is the ^^—diffeomorphism from the 
torus unto dflh- We denote by k the curvature of dO written in local coordinates, 
and let ho be such that 



Asymptotic for an insulating thin layer. The first theorem gives the asymp- 
totic expansion of the solution Vh of |T]), for h tending to zero, for bounded a. 
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Theorem 1. Let h belong to (0,ho). The complex parameter a satisfies 

(3) \a\ < 1, 

(4) $t(a) > or {jt(a) = and 3(a) ^ o|. 

Let N G N and 4> belong to H N+3 ^ 2 (dQh) ■ Denote by f and i the following func- 
tions: 

V6»eM/27rZ, f(6)=4>o<S> 1 (e), 

VxedO, f(x) =4>o^> 1 o% 1 (x). 
Define the sequence of potentials (V k , V k m ) k=0 as follows. We impose 

V(r?,0)eC, W = o, 

and we wse i/ie convention 

U c = 0, ifl<-l, 

\vr = o, ifi<-i. 

For < k < N we define for all < s < 1 i/ie function d v V k m {s, ■) on R/2ttZ : 



«+i(*, o =Ji lfc+ i/ + ^ {« {3^^r + w} 



+ ■ 



V 3 K 3 d 2 v V k m _ 2 + v 2 K 3 d n V£l 2 + n K d 2 9 V k m _ 2 - nK'dgV^dn, 

and the functions V k and V k m are then defined by 
AT4 C - 0, 

d n V k c \ ao = ad T] V k m +l o^\ 
V k c da = 0, 



I dO 

V.s e (0, 1), V k m (s, •) = f d v V k m (n, •) dry + V k c o<S>». 
Jo 

Let R C N and be the functions defined by: 

(R c N = V h -ZL V k c h k ,inO, 
\R% = V h ot>-Y%= V k m h k , in C. 

Then, there exists a constant Cq,n > depending only on the domain O and on 
N such that 

(5a) \\R c n\\hho) < Co,N\\n\H^/2 {d o)Hh N+1/2 , 

(5b) II -Rjv II -Ffi(C) ^ Co,N\\t\\H N +3/2( do ~ ) h N+1/ ' 2 . 

Moreover, if <j> belongs to H N+5 / 2 (dQ.h), then we have 
(6a) \\R C N\\m { o) < CoAnH»+^ (9 o)W\h N+ \ 

(6b) \\Rn\\hI(c) < Co,N\\i\\ H N + 5 / 2 (do) hN+1/2 - 



In this theorem, we approach the potential in the inner domain at the order 
N by solving TV elementary problems with appropriate boundary condition. From 
these results, we may build an approximated boundary condition on dO at any 
order, in order to solve only one problem. However, this kind of conditions lead to 
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numerical unstabilities, this is the reason why we think that the method to obtain 
the potential step by step is more useful. 

Since it is classical to write approximated boundary conditions we make precise 
these conditions at the orders and 1. Denote by & the curvature of dO in Euclidean 
coordinates and by V® pp and V* the approximated potentials with approximated 
boundary condition at the order and 1 respectively. We have: 

(7) 

(8) O.X'L. nf. 

and 



AV° pp = 0, in O, 
d V° 

u n v app 



(9) 
(10) 



AV 1 

r app 

d n vL, 



0, in O, 
■ ahd?V* pp 



a(l + 



where dt denotes the tangential derivative on dO. The boundary condition (fTU|) 
imposed to d„ V^ pp is well-known for non highly contrasted media. It might be 
found in |16j . With our theorem, we prove that it remains valid for a very insu- 
lating membrane, and we give precise norm estimates. Moreover we give complete 
asymptotic expansion of the potential in both domains (the inner domain and the 
thin layer). 

We perform numerial simulations in a circle of radius 1 surrounded by a thin 
layer of thickness h. In Fig'0 the left frame illustrates the asymptotic estimates at 
the orders and 1 of Theorem [T] for an insulating thin layer. However, the right 
frame shows that as soon as the thin layer becomes very conducting, for example 
as soon as a = l/h, these asymptotics are no more valid: we have to use the 
asymptotics of Theorem [5] 

Asymptotics for an insulating inner domain. Let (3 be a complex parameter 
satisfying: 

Re{(3) > 0, or (Re((3) = 0, and 3f(/3) ^ 0) . 
The modulus of (3 may tend to infinity, or to zero but it must satisfy: 



and — - 



" I T" 



h, 

N eN. 



We suppose that a 



Theorem 2. Let h belong to (0,h ). Let q e N* . 
satisfies: 

(11) a = (3h q . 

Let <f> belong to H N+3 ^ 2+q (dQh) o,nd denote by f and f the following functions: 
yeeR/2nZ, f{9) = 0o$i(0), 

yxedO, f(a:)=^o$ 1 o$ 1 (i). 
Define the function (u^' 9 , w™'' 3 )^ = _ 1 by induction as follows, with the convention 

0, ifl< -2, 
= 0, if I < -2. 




Ifq=l 



Au c 



0, in O, 



£)2 c 
Of U_ 



13 d n u_ 



\ uz\&do = 0. 

V J dO 
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(a) a ■ 



10 



10 




o- 

10 -4 

10 



First Order 



10 



10 

* Zeroth Order 



10 



(b) a = i/h 



Figure 2. H 1 norm of the error at the orders and 1. 



Vfa,0)eC, u™? =u c l\\ ao o§ . 



Moreover, 



d v u^ = o, d.uT' 1 = (i - vffiu-i 1 + /• 
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For < k < N, denote by <j>\ the following function: 
4>l = [ (« {^ulX\ + d n u^\) + nndlu^\ nn'dgu^) dn. 



and define u^ 1 by 
' Aui 1 = 0, in O, 



+ P dnuf\ ao = Ui J\r, lJSg^ur 1 ^] o* 1 



dO 

I uf& a o = 0. 

V JdO 

In the membrane u™' 1 is defined by 

u k = / d v u k ' y dr7 + u fe "o$o, 
Jo 

and d^u™^ for i = 1, 2 is determined by: 

n m,l / / /o o2 m,l . o m,l 

• If q > 2. The function u™ ,q is defined by 

[ u m [ q d6 = 0, 

ITie potential u2\ is solution to the following problem: 
f At£? = 0, O, 

Moreover, 

d v u^ = 0, d r ,uT' q = (1 - r?)^ 9 + /• 
For < k < N , denote by (j)f, the following function: 



-l 



- > ! < k (Srjd^u^ll + d v u™±\) + rind 2 e u™l\ - nn'dgu™^ ) drj. 



u™' q \ n =i is entirely determined by the equality: 

Jo 

hence 

The potential u c ^ q satisfies the following boundary value problem: 
f Au£ 9 = 0, in O, 
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The functions (^u^*J ._ satisfies equation (|12p , m which 
placed by u m ' q . 
Let r C pf and r^ ,q be the functions defined by: 

Then, there exists a constant Cq,n > depending only on the domain O and on 
N such that 

\\r c N q \\m{o) ^ Co,jv||f||H« + 3/2 +g(aci) max (J^'^j hN+1/2 , 

II m .9ll ^ r< llrll lJV+1/2 

Fat 1 1 if i(C) S ^0,JV||l||ifJV+3/2( ac ,)ft 

7/0 belongs to H N+5 / 2+ i(dVt h ) , we have 

W r< N q \\m(0) — Co,N\\f\\H N +5/ 2 +i(T)h N+1 ■ 

We observe that if q = 1 and iV = 0, the approximated boundary condition at 
the order is given by: 

-(1 - hft/2)d?ul app + —!^-d t ul app + (3d n ul app = - ^oj'io^ 1 . 

Thus it is very different from the approximated boundary condition |jHJ) imposed 
to V® pp in the case of an insulating membrane. This is a feature of the conducting 
thin layer. Observe on Fig [3] that the numerical computations in a circle confirm 
our theorical results. 




Figure 3. H 1 Norm of the error at the order —1 and for an 
insulating inner domain : a = ih. 



Thanks to our previous results by comparing the parameters \a\ and ft of a 
heterogeneous medium with thin layer, we know a priori, which asymptotic formula 
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(Theorem Q] or Theorem [5]) has to be computed. We emphasize that our method 
might be easily implemented by iterative process as soon as the geometry of the 
domain is precisely known. 

In the following, we show how the potentials (V fe c , V™)k>o of the previous the- 
orems are built, and then we prove these theorems. Let us now make precise the 
geometric conventions. 

1. Geometry 

The boundary of the domain O is assumed to be smooth. The orientation of the 
boundary dO is the trigonometric orientation. To simplify, we suppose that the 
length of dO is equal to 2-k. We denote by T the flat torus: 

T = M/2ttZ. 

Since dO is smooth, we can parameterize it by a function if? of class < ^' 00 from T to 
R 2 satisfying: 

V0eT, |*' (0)| = 1. 

Since the boundary dVlh of the cell is parallel to the boundary dO of the inner 
domain the following identities hold: 

dO = {*(0),0€ T}, 

and 

dn h = {Tf(9) + hn(9),6 eT}. 

Here n(0) is the unitary exterior normal at ^?(9) to dO. Therefore the membrane 
Oh is parameterized by: 

O fe = {$(r?,0), fa, 0) €]0, l[xT}, 

where 

$(r?,0) = ^{6) + hrjn(9). 
Denote by k the curvature of dO. Let ho belong to (0, 1) such that: 

(13) h < — Jj-. 

|| K ||oo 

Thus for all h in [0, ho], there exists an open intervall / containing (0, 1) such that $ 
is a smooth diffcomorphism from / x R/2irZ to its image, which is a neighborhood 
of the membrane. The metric in Oh is: 

(14) h 2 di] 2 + (l + VO 2 d0 2 . 

Thus, we use two systems of coordinates, depending on the domains O and Oh- in 
the interior domain O, we use Euclidean coordinates (x,y) and in the membrane 
Oh, we use local (rj, 0) coordinates with metric (TH]) . 

We translate into the language of differential forms Problem {T]) and Problem 
@. We refer the reader to Dubrovin, Fomenko and Novikov [5] or Flanders [T3j for 
the definition of the exterior derivative denoted by d, the exterior product denoted 
by ext, the interior derivative denoted by 5 and the interior product denoted by 
int. In Appendix we give the formulae describing these operators in the case of a 
general 2D metric. Our aim, while rewriting our problems {1} and |J2J is to take 
into account nicely the change of coordinates in the thin membrane. 

Let V be the 0-form on fih such that, in the Euclidean coordinates (x,y), V is 
equal to V, and let F be the 0-form, which is equal to <f> on dflh- We denote by N 
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the 1-form corresponding to the inward unit normal on the boundary fi/j (see for 
instance Gilkey et al. [H] p. 33): 

N = N x dx + N y dy, 

N* is the inward unit normal 1-form. Problem (JTJ) takes now the intrinsic form: 
(15a) S {q h dV) = 0, in fl h , 

(15b) int(N*)dV = F, on dQ h . 

According to Green's formula (Lemma 1.5.1 of [2]), we obtain the following trans- 
mission conditions for V along dO: 

int(N*)dV| 90 - aint(N*)dV|ao h y?n h , 
(15C) ext(N*)V| ao = ext(N*)V| aOh \ a0fe . 

Similarly, denoting by U the 0-form equal to u in Euclidean coordinates we rewrite 
Problem {2} as follows: 

(16a) S( 7/l dU) = 0, inO A , 

(16b) int(N*)dU = F, on dQ h ; 

the following transmission conditions hold on dO: 

aint(N*)dU| ao = int(N*)dU| 90h \ 9 n h , 
( C) ext(N*)U| S o = cxt(N*)U| a0h \ a n h - 

2. Statement of the problem 

In this section, we write Problem (fT5|) and Problem (fT6|) in local coordinates, 
with the help of differential forms. It is convenient to write: 

V0eT, $ o (0) = $(O,0),$i(0) = $(l,0), 

and to denote by C the cylinder: 

C = [0, 1] x T. 
We denote by / and f the following functions: 

(17) V(x,y)edO, R(x,y) = KO% 1 (x,y), 

(18) V6»gT, f(6)=(t>o$ 1 (6), 

(19) VxedO, f = /o$o 1 (x). 

Using the expressions of the differential operators d and <5, which are respectively 
the exterior and the interior derivatives (see Appendix), applied to the metric (fH| , 
the Laplacian in the membrane is given in the local coordinates (77, 9) by: 

V (77, 6>) S C, 

(20) AL n , e) =T7T ( 1± ^^+ -±-V ' 



I*M) - h ( 1 + hr)K )^ \ h Wr >J^ l + hr]K\l + hr,K 
Moreover, for a 0-form z defined in Oh, we have: 

mt(N*)dz\ go = ^d v z\, 1=0 , 
int(N*)dz| SOh = ^d v z\, 1=1 . 
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Denote by 



V c = V, in O, 
V m = 7o$, in C, 



and by 



u c = u, in 0, 
H m = MO$, in C. 

We infer that Problem (fT5|) may be rewritten as follows: 



(21a) 



AV C = 0, in O, 
1 



(21b) V(77,0)eC, ^((l + ^Kjfl^O + afl 



1 + hrjK 



d e V r 



(21c) 

(21d) 
(21e) 



V c o<S> = V r ' 



r)=0 



\ri=0 ' 



Vda = 0. 



Similarly the couple (u c , u m ) satisfies 



(22a) 



Au c = 0, in O, 



1 



(22b) V ( V , 6) e C, -^d v ( (1 + hr,K)d v u r ' 



1 



1 + /l77« 



0. 



(22c) 

(22d) 
(22e) 



ad n u c o<&o — ~d n U T ' 
h 



r/=0 



u c o$ = u 



\r/=0 ' 

hf, 



i do- = 0. 



Remark 3. In the following, the parameter a is such that: 

5R(a) > or {&(£*) = Oand 9f(a) 7^ o|. 

Since a represents a complex permittivity it may be written (see Balanis and Con- 
stantine [Ijj as follows: 

a = e — ia/w, 

with s, a, and ui positive. Thus this hypothesis is always satisfied for dielectric 
materials. 

Notation 4. We provide C with the metric (|14p . The L 2 norm of a 0-form u in 
C , denoted by ||m||a°l 2 (c)> * s equal to: 

II«IIaol»(c) = ( / / Mi + ^)r"M)l 2 d??d0 



Jo 
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and the L? norm of its exterior derivative du, denoted by ||cIu||ail 2 * s equal to 



|du||AiL2 (C) 



10 JO " 

= l|gradu|| L 2 (0h) . 



1/2 



To simplify our notations, for a 0-form u defined on C, we define by ||w||iri(c?) the 
following quantity 

\\ u \\hI(c) = \\u\\a°lIjc) + l|du|| A i L 2 i(c) , 

when the above integrals are well-defined. Observe that for a function u € H l (Oh), 
we have: 

\\Aw-{o h ) = \\ u0 ^\\hI{C)- 

Remark 5 (Poincare inequality in the thin layer). Let z belong to i/r (C), such 
that 



(23) 



/■27T 

/ z(O,0)d0 = O. 
Jo 



Then, there exists an h-independant constant Co such that 
( 24 ) IM|a°£3(c) ^ °o l|dz|| A i L 2 (c) ■ 

We prove (|24p using Fourier analysis. According to the definition (| 13[) of ho there 
exists two constants Co and cq depending on the domain O such that the following 
inequalities hold: 

-1 /.27T 



(25a) \\z\\io LHD) < C Q h / / \z<ji,0yC M&T) 



o Jo 



(25b) 



M\vlI{d) > c o 



\d v z(r),i 



h\d z\ 2 d6dn 



For k G Z, we denote by the k th -Fourier coefficient (with respect to 9) of z: 



1 r 2 * 

Zk = T- z(9)e- ike d9. 
2ti" Jo 



Since [dgz] — ikzk, we infer: 



Vfc^O, / \z k (v)rdrj< 
Jo 

According to gauge condition (|23p . we have: 

z (0) = 0, 
</ims ; using the equality 



(**) fo) 



drj. 



we m/er 



|2b(?y)l dry < / (<VJ (»?) 



dr/. 
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Therefore, 



£ /V*M)l a d*7<£{ f\(?**) M 



drj 



(v) fc fo) 



We end the proof of 



by using Parseval inequality and inequalities ([2"5 



3. Asymptotic expansion of the steady state potential for an 

insulating membrane 

We derive asymptotic expansions with respect to h of the potentials (V c ,V m ) 
solution to Problem (|2ip . The membrane is insulating since the modulus of a is 
supposed to be smaller than 1. However, our results are still valid if \a\ is bounded 
by a constant Co greater than 1. We emphasize that the following results are valid 
for a tending to zero. 

3.1. Formal asymptotic expansion. We write the following ansatz: 
(26a) V c = V Q C + hVf + h 2 V 2 c + ■ ■ ■ , 

(26b) v m = v Q m + hvi n + h 2 v; n + ■■■ . 

We multiply (|21b|) by h 2 (l + hnn) 2 and we order the powers of h to obtain: 

V(?y, 9) e [0, 1] x T, 

d 2 V rn + hn { 3r]d 2 V m + d v V m ) + h 2 {3n 2 K 2 d 2 V m + 2nK 2 d v V m + d 2 e V rn } 

' " ' + h 3 {r} 3 K 3 d*V m + n 2 K 3 d v V m + nnd 2 e V m - vk 1 ' d e V m ) = 

We are now ready to derive formally the terms of the asymptotic expansions of V c 
and V m by identifying the terms of the same power in h. 

Recall that for (m, n) in N 2 , 6 m ,n is Kronecker symbol equal to 1 if to = n and 
to if to 7^ n. By identifying the powers of h, we infer that for I £ N, Vf and V™ 
satisfy the following equations: 

(28a) AVf =0, in O, 

for all (r?, 6) € C, 



d 2 vr 



(28b) 



+ 3n 2 K 2 d 2 V^ 2 + 2n K 2 d v V^ 2 + d 2 Vfl 



V,': 



with transmission conditions 

(28c) d n V l c o^ =ad v V l X 1 \ ri=0 : 
(28d) V l c o$ = Vr\ n = , 

with boundary condition 
(28e) d v Vr\ n=1 =S h if, 
and with gauge condition 



(28f) 



/ Vfda =0. 

JdO 
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In equations (j28|) . we have implicitly imposed 

r = 0, if I < -1, 

(29) ( v; n = 0, if I < -1. 

The next lemma ensures that for each non null integer N, the functions V$ and 
V™ are entirely determined if the boundary condition (f> is enough regular. 

Notation 6. For s G R, we denote by < ^' 00 ([0, 1]; H S (T)) the space of functions 
u defined for (rj,6) G [0,1] x T, such that for almost all 8 G T, u(-,6) belongs to 
■af 00 ([0,1]), and such that for aline [0,1], u(r),-) belongs to H S (J). 

Lemma 7. We suppose that DO is smooth. 

For N G N and p > we suppose that cf> belongs to H N+p ~ 1 / 2 (dQh) and let 
\a\<l. 

Then the functions V m , • • • , Vff and V§, ■ ■ ■ , V§ are uniquely determined and 
they belong to the respective functional spaces: 

Vfe = 0, • • • , N, 

(30a) Vl n G ([0, l];H N+p - k+1 / 2 (T)) , 

(30b) V£ G H N+p - k+1 (0). 

Moreover, there exists a constant Cn,o.p such that: 

Mk = 0, • • • , N, 

(31a) sup 11^(77, OH 

ne[o,i] 

(31b) \\V k c \\ n N+p-fc+i(Q) < Ic^lCv.cpllf ll_f/ N +p- i /2(ac>)- 



Remark 8. To simplify, we suppose that \a\ < 1, but the same result may be 
obtained if there exists Cq > 1 such that \a\ < Cq. In this case, the constant Cn,o,p 
would also depends on Cq. 

Proof. Since dO is smooth and since 4> belongs to H N+ P~ 1 / 2 (dn h ), for N > and 
p > 0, then the functions / and f defined by (fl8]) and by (|19| belong respectively to 
i? Ar+p_1 / 2 (T) and to H N+p ~ 1 / 2 (dO). We prove this lemma by recursive process. 
• N = 0. Let p > and let (j) belong to H p -^ 2 {dQ. h ). 

Thus / and f belong respectively to H^ 1 ' 2 ^) and Ip-^dO). Using (l28bl) and 
Sell , we infer: 



, WU = o, 



(32) 

hence, c^Vq" 1 = 0. According to (|28b[) and to (|28c[) . we straight infer 

w = /. 

Therefore by (|28a[) and (|28c|l the function V C satisfies the following Laplace prob- 
lem: 

(33a) AV C = 0, 

(33b) d n V c \ dO = at, 

with gauge condition 

(33c) / Vo C dcr = 0. 
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According to (|28d|) . we infer 
(34) ^ m = ^ c o$ , 

hence Vq and Vg™ are entirely determined and they belong to the following spaces: 

F m e ([0,l];iP +1 / 2 (T)) , 
V£ e H p+1 (0). 

Observe also that there exists a constant Co, P such that 

sup \\Vq 1 (t] 1 -)\\ Hp+1/ 2 (t) < Co, P \\{\\HP-wtao)) 

\\Vo\\hp+i(0) ^ \ a \Co,p\\{\\HP-i/2(dO)- 

• Induction. 

Let N > 0. Suppose that for all p > 0, for all <f> E H N+ P- 1 / 2 {dfl h ) and for 
M = 0, • • • , N the functions V£ f and V]$ are known. Suppose that they belong 
respectively to H n+ p- m+1 (0) and to £ c g°° ([0, 1]; H N+p - M+1 / 2 (T)) and that 
estimates (j3"Tj) hold. 

Let belong to H N+v+1 ' 2 (dQ. h ). Therefore, for M = 0, ■ • • , N, the functions 
Vfa and V$ are known, they belong respectively to h n +p- m+2 {0) and to V$ € 
c g°° ([0, 1];H n +p- m+3 / 2 (T)) and the following estimates hold: 

VM = 0, • • • , N, 
SU P \\Vm (v> OlljjJv+p-M+s/arn — ^,o,p||f||ffw+ P +i/2( a0 ), 

»)G[0,1] 

II^Mll_f/«+p-M+2(ci) — l a |CA r ,Cpl|f|lif JV +P+ 1 /2(9c>)- 

We are going to build V^ +1 and V^ +1 . From (|28b|) and (|28e[) . we infer, for all 

a?v^+i = - 1« {3^yxp + w} 

+ r, 3 « 3 9 2 y™_ 2 + t? 2 * 3 W_ 2 + rfl^_ 2 - ^flsV^aj, 



Recall that we use convention (|29|) . Since we have supposed that VjJJ? is known for 
M < N and belongs to ([0, 1]; h n+1+ p- m - 1 ' 2 (T)) , we infer that: 

V(M)eC, 

(35) 

W+i(*> •) =J {« {3»^V^ + W} 



ASYMPTOTICS FOR STEADY STATE VOLTAGE POTENTIALS 



17 



is entirely determined and belongs to < t?°° ([0, 1]; ff p+1 / 2 (T)) . Moreover, since 
d 71 V^ l +1 is known, we infer exactly by the same way that d^V^^ is also deter- 
mined. Actually, it is equal to 

V( s ,0)ec, 

+ fjVflgV^! + ti 2 k%VS>L 1 + ^IVn-i ~ ^Wjv-i f dry, 



and it belongs to < ^°° ([0, 1]; iJ p+1 / 2 (T)). According to pgcj) . the function V£ + i is 
then uniquely determined by 



(36a) AV£ +1 = 0, 

(36b) W +1 | ao = aft,V?P +2 o*o S 

with gauge condition 

(36c) / V$ +1 6tr = 0. 

JdO 

Moreover, it belongs to H p+1 (0). Transmission condition (|28d| implies the fol- 
lowing expression of V^ +1 : 

Ms e (0, 1), V$ +1 (8, •) = f W+ifa> •) drj + V£ +1 o$ , 

Jo 

where ^V^ +1 is given by (35]) and belongs to ([0, 1]; HP+^ 2 (T)) . We infer 
also that there exists Cn+i,o.p > such that 

sup || V^ l +1 (r], ■)\\ HP+1/2m < CV + i iC i,p||i||ff« +P+ i/2( 0O ), 
»?e[o,i] 

IIVSv+iIIhs+^O) - \ 0l \ < ~'N+l,0,p\\t\\H N +P+ 1 / 2 {dO)i 

hence the lemma. □ 

Observe that the functions (V fc c , V™) are these given in Theorem [TJ 

3.2. Error Estimates of Theorem [TJ Let us prove now the estimates of Theo- 
rem UJ Let N e N and <f> belong to H N+3 / 2 (dn h ). The function f is defined by 
(fT5| . Let and R™ be the functions defined by: 

\R% = V h o$-j:Z =0 V k m h k ,m C. 

We have to prove that there exists a constant Cq,n > depending only on the 
domain O and on N such that 

(37a) \\R c n\\huo) < C ,N\\n\H»+^(do)\a\h N+1/2 , 

(37b) \\Rn\\hi(c) < Co,v||i||H«+3/2( ac i-)/i Ar+1//2 . 

Moreover, if <j> belongs to H N+5 / 2 (dn h ), then we have 
(38a) WRnWhHO) < Co,N\n H "^ /Hdo) \a\h N+ \ 

(38b) \\R"n\\hi(c) < Co,v||f||H« + 5/2( 9c ,)/i Ar+1/2 . 
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Proof of Theorem\7\ Since <ft belongs to H N+3 ^ 2 (dilh), according to the previous 
lemma, the couples of functions {R%, and (i2^ +1 , RJf +1 ) are well defined and 
belong to H l {0) x H\ (C). The Sobolev space H] (C) is defined in Notation H 
Denote by g^ the following function defined on C: 

g N =k (3r)d*V£ + d v Vtf) + 2,if k 2 d 2 ^^ + 2r l K 2 d r ,Vtf_ 1 + d 2 g V^_ x 

+ ^^d 2 v V^_ 2 + V 2 K 3 d v V^_ 2 + wd 2 V^_ 2 - r,K'd e V^_ 2 



(39) 



+ h \^ v 2 K 2 d 2 v^ + 2 V K 2 d v v£ l + d 2 g vi? 



+ h 2 rf^d 2 ^ + v 2 k%V^ + T]Kd 2 B V^ - ■ n K , d a V 1 



N 



According to the previous lemma and since <fi belongs to H N+1 ' 2 (d£lh), the above 
function g N belongs to ([0, 1]; ii" 1/2 (T)) and the function fyVff belongs to 
([0, l]; J ff 3/2 (T)). Moreover, there exists a constant C N .o such that 



(40) 



SU Pr)e[0 : l] ll.9A r (^; -)llff-i/2(T) < C N ,0 \\ f \\ H « + V2(T) ) 
v SU Pr,G[0 : l] \\ 9 r, V N (Vi -)llif3/2( T ) < CjV,C>||/|Un+1/2(T)' 

The functions R C N and R^ satisfy the following problem: 
AR C N = 0, in O, 

with transmission conditions: 

d n R c N o$ = \ (d v R%\ n=Q + h N+1 d v Vj 

R c N o®o= Rn\t,=o > 
with boundary condition 

d V R N\ v= i = 0, 

and with gauge condition 
R C N da = 0. 



n v n \ v =o I i 



By multiplying the above equality by Rn and by integration by parts, we infer 
that: 



c 



(41) +*h N+1 I d v V^\ v=0 R%\, 1=0 d9 



T 



By hypothesis ([4]), and using by Cauchy-Schwarz inequality and estimates ()40|) . we 
infer that there exists a constant Cq,n > such that 

Sft(a) ||di2jvl| A i£3 (C7) - l a l G 0,A r ^ _1/2 |I^JvllffJ(C')ll/ll_H-« + i/2( T )||i?^||_ H 'i(c), 
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and 

|9f(a)| \\dR%f AlLl{c) < lalCc.iv^-^ll^lk^oll/ll^+v^Dll^llH^c), 
hence 

\\^ R n\\mlI(C) - C 0,Nh N ~ 1/2 \\ f\\ h n + i/2(j) || RnWhi(c) ■ 

Since / T Rf}\ v=0 d9 = 0, by Poincare inequality (|24|. there exists a strictly positive 
constant Co, which does not depend on h such that 

\\Rn\\a°LI(C) < Co l|di?^|| A i L 2( C) , 

hence 

\\ r n\\hi(c) ^ Co,v||/||ff« + i/2( T) /i Ar_1/2 , 
and therefore we deduce directly from the above estimate and from (j¥T|) , 
\\ r n\\h^{o) ^ Co,v||/||ii-« + i/2 (T) |a|/i Ar_1/2 . 

The above estimate holds for E H N+1 / 2 (dfl h ). Since belongs to H N+3 / 2 (dQ h ), 
we obtain the same result by replacing N by N + 1: 

|-Rjv+i||ffi(c) — Co,v+i||/||HN+3/2( T )/i Ar+1 / 2 , 

11^1(0) - C 0,N+l\\f\\ H N + :i / 2 (T)\ a \ hN+1/2 - 

According to the previous lemma, the functions Vfc +1 and V^p +1 are well-defined 
and there exists a constant Cm,o such that: 

\\ v n+i\\h^(o) < |"|Civ,ei||f 11^+3/2(90), 

II^JV+lll-ffJ(C) < -^-||f||//«+3/2(9C))- 

^iV — + ^V+1^ JV+1 ' 

nm pm , t/to i.JV+1 

/% — JtjV+1 "+" • / V+l n ' > 

(T)l«|^ +1/2 , 



(42) 



Writing 
and 

we infer that 
and 



II^Hfl-jCC) - Co,Ar||/|| jH -«+3/2( T )/i 7V+1/2 . 
If (f> belongs to H N+5 ^(dQ h ), we write 

7->e _ 7->c | ire h N+1 4- T/ c ?7 Ar + 2 

"JV — JX N+2 T ^Ar+l' 4 T • / AT+2' t ) 

and 

n N — n N+l T t'jV+1' 1 ' ,/ V+2"- ) 

to obtain estimates hence Theorem [TJ □ 
Remark 9 (The case of an insulating inner domain). Consider Problem 

div (7/j gradu) = in fi^, 

^ = 4> on dfl h , 
on 

u da = 0. 

so 
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// the inner domain is perfectly insulating (i.e. if jh vanishes in O), the steady 
state potential in the membrane satisfies: 



= 0, in C, 



((1 + hr ]K )d v u m ) + d e ^ ^L— frt 
with the following boundary conditions: 

^« m Uo=°. d^"\=i = hf, 
and with gauge condition 

2tt 

u m \ v=0 d6 = 0. 

By identifying the terms of the same power of h we would obtain: u™ = 0, and u™ 
would satisfy: 

d^u™ = 0, in C, 
«l, )= o = 0, d v u?\ v=1 = f, 

u[ n \ v=0 d9 = 0, 



which is a non sense as soon as f 0. Our ansatz (|26p fails. Actually, the 
asymptotic expansion of u m begins at the order — 1 : a boundary layer phenomenon 
appears. This is described in the next section. 

4. Asymptotic expansion of the steady state potential for an 
insulating inner domain 

Consider now the solution Uh to Problem {2}. In this section, we suppose 

(43a) \a\ tends to zero, 

(43b) 5R(a) > or jsR(or) = and 3(a) ^ oj. 

Thus the inner domain is insulating. Let f3 be a complex parameter satisfying: 

Re(f3) > 0, or (Re(/3) = 0, and S(/3) 7^ 0) . 
The modulus of j3 may tend to infinity, or to zero but it must satisfy: 

— I , and - — - = o ( — 

hj \/3\ \h 

We suppose that u may be written as follows: 

Uh = t u -i + u o + hui + ■ ■ ■ . 
h 

We denote by u c and w m o$ _1 the respective restrictions of Uh to O and to Oh- 
One of the two following cases holds. 

Hypothesis 10 (a — (3h q ). There exists q > 1 such that: 

(44) a = (3h q . 

Hypothesis 11 (a = o(h N ), VW € N). TTie complex parameter a satisfies (|43|) 
and /or a// TV S N. 

(45) ViV e N, |a| = o(^). 
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First we suppose that Hypothesis IT01 holds : we will discuss on Hypothesis [TTIlater 
on. We denote by (u°' q ,u m ' q ) the solution to Problem [5] under the Hypothesis fTtJl 

According to |22|) . by ordering and identifying the terms of the same power of 
h, for k £ NU -1, for q € N*, <' 9 and u™' q satisfy: 

(46a) Au^ q =0, in O, 

for all (r?, 6) £ C, 

<9 2 u™' 9 = - |k {3r)d*u^!l + d v u T-i} 
(46b) + 3t7 2 k 2 92 u ™.| + 2?7K 2 3 r) M™'| + d 2 e u^ q 

(46c) «?' 9 o$o = «rVo> 

(46d) W\=i =*,,!/, 

(46e) / u^da =0. 

Transmission condition (|22c[) coupled with Hypothesis [TUl implies: 

(46f) R^ 1 _ g o$ = W\=o. 
In equations (fH)]) . we have implicitly imposed 

jw^ 9 = 0, if I < -2, 
|up 9 = 0, if/<-2. 

Let us now derive formal asymptotics of u when Hypothesis 1101 holds. 



(47) 



4.1. Formal asymptotics. 

• N = -1. 

The functions satisfies 

f 9 2 U m i 9 = 0,in C, 

\ ^u™i Q, |^=o = 0, 9^n_ T '| 9 |^ = i = 0, 

hence u™[ q depends only on the variable 0. Observe that we have, for almost all 
6 6 T the following equality: 

• N = 0. 

The function it™' 9 satisfies: 

f3X' 9 = 0,in c, 

[ 9^u^' ? |^=o = 0, d n u™' q \ n =\ = 0, 

hence, d v u™' q vanishes identically in C. 

• N = 1. 

The functions w™' 9 satisfy: 

J ^-u™' 9 = -<^u™i 9 ,m C, 

[ 9^u7' 9 |,,=o = PdnU'zl _ q o $0, i9 J) w™ ,9 |^ = i = /. 
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Therefore for q — 1 we obtain the following equality: 

-dgu^ 1 + (3d n uz\o<& a = f, 
hence the following boundary condition imposed to uz\ on dO: 

= f. 



£)2 c,l 

d ul t 



DO 



+ {3 d n u c l\ 



DO 



Therefore, the function u c j\ is solution to the following problem: 



' Au c _l\ = 0, in O, 



(48) 



(49) 



■~i2 c,l 



oo 



P d n u c l\ 



oo 



/ u c l\& do 

JdO 



u c l\\doo§ . 



Since 5R(/3) > 0, a straight application of Lax-Milgram theorem ensures that zi Cl1 is 
uniquely determined and belongs to H l {0) as soon as the boundary data belongs 
to H-*/ 2 {dO). 

If q > 2, the function u™' 9 satisfies: 



(50) 



2 u m \ q = 



Since JjU^fdO = 0, equality (|50| defines uniquely We infer that uz\ is 

solution to the following problem: 



(51) 



,c,q _ 



0, in O, 



Hence we have determined it™{ 9 and u^f for q € N*. Observe that uZ^ is solution 
to Laplace equation with mixed boundary condition, and for q > 2 the potential 
u^zf is the solution to Laplace equation with Dirichlet boundary condition, while for 
an insulating membrane, we obtained Neumann conditions for the approximated 
steady state potentials. 



Let us now determined 
• Induction. 



m,q 
1 N 



and Upf for q £ N* by recurrence. 



Suppose that for N > 0, the functions u N ,q 1 



The function 



L N+2 



•*7V — 1' 9 V U N 



satisfies: 



V rn, q 



»2 m,q 

OgU N 



d2 m,q i W,Q • /^t 

, e u Ar_i + ?7 K Oeu N _ 1 , m C, 



Denote by <j> q N the following function: 



k (3rjd 2 



N+l 



u V a N+l) 



,2 oi m,ij 
N-l 



7]K'dgU™[ q j d.77 



Since d^u 7 ^'^ and d^u™'^ are supposed to be known, the function <fi q N is entirely 



determined. Observe that if q — 1, 9^w^^ 2 | ??= o is unknown since d n u c ^ is not yet 

m,q I 
JV+2 l'7 ; 

<•) , we infer the following equality satisfied by 



determined, while as soon as q > 2, d v u™'£ 2 \r]=o is known. 



Using transmission condition 

m.l • n 

u N m r\ = 0: 



d2 m,l I 
\V=0 



PdnV.% 1 O $0 =0]y - / 
JO 
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hence the boundary condition imposed to u c ^ on dO 



P d n u 



N 



DO 



£)2 c,l 



Thus the function u^ 1 is solution to the following problem: 



r Au 



(52) 



V 



0, in O, 



a2 c,l 



ye> 



13 d n u 



N 



i)0 



4>N- / (V-I^d^dn 0$o 



u^dao = 0. 



In the membrane u^' is defined by 

(53) "tv' 1 = / d v u^' g dri + u% q o^o- 

Jo 



If <7 > 2, Ujy 1,7=1 is entirely determined by the equality: 



-dju^\ =1 =/3d n u% q + x _ o$ + </>%-! (V- ^d'edrju^dr), 



hence 



The potential u^ 9 satisfies the following boundary value problem: 



(54) 



A' 



c,g _ 



0, in O, 



% lyO = U N 0$ n 



Observe that for g > 1, d n u N '? 2 is then entirely determined by: 



K (377^ 



dew^' 9 - V Kd e u N-i + riddeu™'^ dry. 



Therefore, we have proved that for all N > —1, for q G N*, the functions and 
w^' 9 are uniquely determined. 

Remark 12 (Regularity). Observe that these functions are the potentials given in 
Theorem [Jl We leave the reader verify by induction that the following regularities 
hold. Let q € W, N > -1 and p > 1. Let <j> belong to H N+ P- 3 / 2 (dflh). 

u c l\ e H 1+N+ '" (O), 

u™{ q e ([0, l];fl-V3+^+P(T)) , 



VA; = 0, • • • , N, 



(55a) 
(55b) 



■ q G < ^ >0 ° (j0, l];# 1/2+Ar+p ~ 2[(fe+1)/2l (T 
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Moreover, there exists a constant Cn,o,p independant on h and (3 such that: 



(56a) sup •)|| ff i/2 + N +PfT) < Civ,o,p||f||ifJv+p-3/2(aci), 

»?e[o,i] " y ' 

(56b) < C < A',o,p||i||ifN+p-3/2( ac i), 

Vfc = 0, • • • , N, 

(56c) sup \\u™' q (r], -)llfl-i/2+iv+p-2[(fe+i)/2] m < Cjv,o,p||f||ffJv+p-3/2(90), 

(56d) ll M fc 9 llHi+Jv+ P -[fc/2](c)) - CA r ,o,pPllff JV +p- 3 / 2 (ac')- 



4.2. Error estimates of Theorem [2} Let us now prove Theorem[3J Let q E N* 
and JV 6 N. The complex parameter a satisfies (|4"5|) with Hypothesis [TUJ Let 
belong to H N+3 / 2+q (dflh)- Let r^ 9 and r^' 9 be the functions defined by: 

We have to prove that there exists a constant Ce>,jv > depending only on the 
domain O and on TV such that 

(57a) \\r% q \\m(o) < CoM^Wh^+^+i^o) max 

(57b) |k™' 9 ||^ (C) < C a Ar||f|| H « + 3/2 (ao) ^ +1 / 2 . 

If </> belongs to ff JV+5 / 2 +9(af2/ l ) ) we have 

H r ;v 9 llffi(0) - Co,Jv||/||ffN+5/2+,( T) /i Ar+1 . 

Proof. The proof of Theorem [5] is similar to the proof of Theorem [TJ Since </> 
belongs to H N ~ 1 ' 2 (d£lh), according to the previous lemma, the couples of functions 
(r^ 9 ,r^' 9 ) and (r^f + 1; r^i) are well defined and belong to iJ 1 (C) x if* (C). 
Denote by <?j\r the following function defined on C: 

g N =k (3rjd 2 u^ q + d n u^ q ) + 3rj 2 k 2 dfa™'^ + 2r]K 2 d v u^' q 1 + d^u^'li 
+ r) 3 K 3 d 2 u^' q 2 + rj 2 K 3 'd v U r ^ [ % + r\Kd 2 u^' q 2 - tjk 1 ' d s u^' q 2 



+ h (^ V 2 K 2 d 2 u^ q + 2r, K 2 d v u™' q + d 2 u^ q 
+ i] 3 K 3 d 2 u^- q 1 + v 2 K 3 d ri u^ q 1 + r]KdgU^ q 1 - r}K,'d e u™' q ^J 

+ h 2 ^r] 3 k 3 d 2 u r £' q + rj 2 n 3 d ri u^ q + r]Kd 2 u^' q ~ rjn'deu^ 1 ^ 

According to the previous lemma and since <fi belongs to ff Ar_1 / 2 (9n? ! ,), the above 
function belongs to < ^°° ([0, 1]; iJ _1 / 2 (T)) and the function d v V^ l ' q belongs to 
< ^°° ([0, 1]; iJ 3 / 2 (T)) . Moreover, there exists a constant Cn.o such that 



(58) 



(59) 



| su P77e[04] \\9n(v, Olln-vafT) < CAf,o||/||irJv-i/2 (T) , 
|sup^ e[0jl] \\dr ! u 7 ^' q {ri,-)\\ H 3/2 m < Cjv,£)||/||h^-V2(t)- 
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The functions r c ^ q and r^' q satisfy the following problem: 
Ar% q = 0, in O, 

1 + hrjK a ^ m ^ q \ , Q ( h a _ miq \ -h N 



with transmission conditions: 
0h 1+ «d n r c jf o$ = ~ (d v r™X =Q + f3h N+1+q fan^U o $ + hd v u™' q o $ )) , 

I M ° w — 'A 1,7=0 ' 

with boundary condition 

Ur l N l)j=l u ' 

and with gauge condition 

By multiplying the above equality by rjy and by integration by parts, we infer that: 
(60) 

Ph'+ q \\dr%% lL2(0) + ||dr">*||^ £S(c) = -h N [j N ( V ,e)r^(r),e)d V de 

d0. 



+ (3h N+1+q / (0„Utf-i o $ + hd v u% q o $ ) r™' 1 

The end of the proof is similar to Theorem [TJ Using the positivity of we 
straight infer estimate (|57b|) of r^' q . To obtain the estimates of r^ q , we write: 

1 

c, Q _ c,g , \ "* c,q ,N+k 
'N — 1 N+q^ 2-^1 N+k 
k=l 

□ 

4.3. The case a = o(h N ), ViV 6 N. Now, we suppose that Hypothesis [TT1 holds. 
In this case, we prove that u c and u m may be approximated by U c and [/ m , which 
are solution to: 

(61a) AU m = 0, in O h , 

(61b) d n U m \ ao = 0, ^f/" l | 9f2h = 0, 

(61c) / C/ m dcr = 0. 

JdO 

and 

(62a) AU C = 0, in O, 

(62b) U% =U m \ 90 . 
Actually, we have the following lemma: 

Lemma 13. Let <j> belong to H^ 1 / 2 (dVlh). Let (u c ,u m ) be the solution to Prob- 
lem and U m and U c be defined respectively by (|61[) and (|62[) . Then, we have: 

(63) \\u m - U m \\ Hl[0h) < C \a\ \<P\ H -y 2{dnh) , 

(64) |K - U c \\ m{0) < CoVW\\<t>\ H -^ { m h ) ■ 
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Proof. Denote by w c and w m the following functions: 

w c = u c -U c , w m ^u m -U m 1 
and let 4> belong to H~ 1 / 2 {d^lh)- We have: 
(65a) Aw c = 0, in O, 

(65b) Aw"' = 0, in O h , 

(65c) a d n w c \ do = d n w rn \ do - a d n U c \ ao , 

(65d) ^\ go = w m \ ao , 

(65e) d v w m \ dnh =0, 

(65f) / w m dcr = 0. 

JdO 

Thus we infer: 

(66) a I \Vw c \ 2 dvo\ + [ \Vw m \ 2 dvo\ 0h = a [ d n U c \ do lF^da. 

JO JO h JdO 

It is well-known that : 

W um \\m(o h ) ^ c o klff-i/2(ao h ) > 

and 

ll^ c || ff i(O h ) < Co \U m \ do \ H i/2 {do) ■ 
Since a satisfies PS| we infer, 

\\w m \\m(o h ) ^ °oH \4>\ B -i/*(on h ) » 

and thereby 

\\w c \\m(o) ^ c oVW\ \<t>\ H -v*( d n h ) ■ 

□ 

It remains to derive asymptotics of U m and then these of U c . They are similar 
to asymptotics of u m ' q for q > 2: we just have to replace (5 by zero. We think the 
reader may easily derive these asymptotics from our previous results. 

Conclusion 

In this paper, we have studied the steady state potentials in a highly contrasted 
domain with thin layer when Neumann boundary condition is imposed on the ex- 
terior boundary. We derived rigorous asymptotics with respect to the thickness of 
the potentials in each domain and we gave error estimate in terms of appropriate 
Sobolev norm of the boundary data, electromagnetic parameters of our domain and 
a constant depending only on the geometry of the domain. It has to be mcntionned 
that for an insulating inner domain (or equivalently a conducting membrane), the 
asymptotic expansions start at the order -1 and mixed or Dirichlet boundary con- 
ditions has to be imposed on the asymptotic terms of the inner domain. 

To illustrate these asymptotics, numerical simulations using FEM are forthcomig 
work with Patrick Dular from Universite de Liege and Ronan Perrussel from Ampere 
laboratory of Lyon. Few results have been shown at the conference NUMELEC 
[TS] with GetDP[TD]. The main difficulty in illustrating the convergence of our 
asymptotic consists in the geometrical approximation of the domain: high-order 
geometric elements seem to be necessary. 
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Appendix 

Let * denote the Hodge star operator, which maps 0-forms to 2-forms, 1-forms 
to 1-forms and 2-forms to 0-forms (see Flanders [E]). We give explicit formulae for 
the operators d, <5, ext and int. These formulae are straightforward consequences 
of the definition of the operators *, d and S — * -1 d*. We refer the reader to 
Dubrovin, Fomenko and Novikov [5]. 

We consider the metric given by the following matrix G 



(67) G 



9n 9i2 

912 922 



We denote by \G\ the determinant of G. The inverse of G is denoted by G 1 

G- 1 = (g^h, 

and we suppose that the signature of G is equal to 1. Thereby, the operator * 2 is 
equal to Id on the space of 0-forms and 2-forms and it is equal to — Id on 1-forms. 

4.4. Star operator in R 2 . 

4.4.1. On 0-forms and on 2-forms. Let T be a 0-form and let S be the 2-form 
vdy 1 dy 2 . Then -kT is the 2-form /idy 1 dy 2 and *S is the 0-form /. The following 
identities hold: 



f= 1 



4.4.2. On 1-forms. Let T be the 1-forrn Tidy 1 + T 2 dy 2 . Then *T is the 1-form 
Hi dy 1 + [i2 dy 2 , and we have the following formulae: 



^ = -V\G\(g 12 T 1+ g 22 T 2 ) 
m = VW\(9 11 Ti+g 12 T 2 ). 



4.5. The action of d acting on 0-forms in R 2 . Let /ibeaO form, then d/i has 
the following expression: 

9/x x dfl 2 
d/i = -z-rdy + -^dy . 
oy 1 oy z 

4.6. The action of S acting on 1-forms on M 2 . Let /i be the 1-form /iidy 1 + 
fi 2 dy 2 , and define 5[i — a. The 0-form a is equal to: 

1 r 9 -^/^(. 9 i Vi+.9 i v 2 ; 



\G\\d yi \ 

+ ^(^(^ + ^ 2 ))}- 



4.7. The exterior product of a 1-form with a 0-form. Let N be the 1-form 
iVidy 1 + N 2 dy 2 and / be a 0-form. The exterior product of ext(JV)/ is: 

ext(iV)/ = fN^y 1 + fN 2 dy 2 . 

4.8. The interior product of a 1-form with a 1-form. Let N and fj, be the 

1-forms iVidy 1 + N 2 dy 2 , and /iidy 1 + fi 2 dy 2 . Then 0-form mt(N)/j has the following 
expression: 

mk(N)ft = m (mg 11 + M2.9 12 ) + ^2 (^ig 12 + V2g 22 ) ■ 



28 



CLAIR POIGNARD 



References 

[1] Habib Ammari and Sailing He. Effective impedance boundary conditions for an inhomogc- 
neous thin layer on a curved metallic surface. IEEE Transactions on Antennas and Propa- 
gation, 46(5):710-715, 1998. 

[2] Habib Ammari and Hyeonbae Kang. Properties of the generalized polarization tensors. Mul- 
tiscale Model. Simul., l(2):335-348 (electronic), 2003. 

[3] Habib Ammari and Hyeonbae Kang. Reconstruction of conductivity inhomogeneities of small 
diameter via boundary measurements. In Inverse problems and spectral theory, volume 348 
of Contemp. Math., pages 23-32. Amer. Math. Soc., Providence, RI, 2004. 

[4] Balanis and Constantine. Advanced Engineering Electromagnetics. John Wiley and Sons Ltd, 
1989. 

[5] Elena Bcretta, Elisa Francini, and Michael S. Vogelius. Asymptotic formulas for steady state 
voltage potentials in the presence of thin inhomogeneities. A rigorous error analysis. J. Math. 
Pures Appl. (9), 82(10):1277-1301, 2003. 

[6] Elena Beretta, Arup Mukherjee, and Michael Vogelius. Asymptotic formulas for steady state 
voltage potentials in the presence of conductivity imperfections of small area. Z. Angew. 
Math. Phys., 52(4):543-572, 2001. 

[7] Yves Capdeboscq and Michael S. Vogelius. A general representation formula for boundary 
voltage perturbations caused by internal conductivity inhomogeneities of low volume fraction. 
M2AN Math. Model. Numer. Anal., 37(1): 159-173, 2003. 

[8] B. Doubrovine, S. Novikov, and A. Fomenko. Geometrie contemporaine. Methodes et appli- 
cations. I. "Mir", Moscow, 1982. Geometrie des surfaces, des groupes de transformations et 
des champs. [Geometry of surfaces, groups of transformations and fields], Translated from 
the Russian by Vladimir Kotliar. 

[9] B. A. Dubrovin, A. T. Fomenko, and S. P. Novikov. Modern geometry — methods and appli- 
cations. Part I, volume 93 of Graduate Texts in Mathematics. Springer- Vcrlag, New York, 
second edition, 1992. The geometry of surfaces, transformation groups, and fields, Translated 
from the Russian by Robert G. Burns. 
[10] P. Dular, C. Geuzainc, F. Hcnrotte, and W. Lcgros. A general environment for the treat- 
ment of discrete problems and its application to the finite clement method. 34(5):3395-3398, 
September 1998. 

[11] B. Engquist and J.C. Nedelec. Effective boundary condition for acoustic and electromagnetic 

scattering thin layer. Technical Report of CMAP, 278, 1993. 
[12] E.C. Fear and M.A. Stuchly. Modeling assemblies of biological cells exposed to electric fields. 

IEEE Trans. Bio. Eng, 45(1):1259-1271 (electronic), 1998. 
[13] Harlcy Flanders. Differential forms wih applications to the physical sciences. Academic Press, 

New York, 1963. 

[14] Peter B. Gilkey, John V. Leahy, and Jconghycong Park. Spectral geometry, Riemannian sub- 
mersions, and the Gromov-Lawson conjecture. Studies in Advanced Mathematics. Chapman 
& Hall/CRC, Boca Raton, FL, 1999. 

[15] K. Idemcn. Straightforward derivation of boundary conditions on sheet simulating an 
anisotropic thin layer. Electron. Lett, 24(ll):663-665, May 1988. 

[16] L. Krahenbiihl and D. Mullcr. Thin layers in electrical engineering, example of shell mod- 
els in analysing eddy-currents by boundary and finite element methods. IEEE Trans. Mag., 
29(2):1450-1455, 1993. 

[17] C. Poignard. Methodes asymptotiques pour le calcul de champs clcctromagnctiqucs dans des 
milieux a couche mince, application aux cellules biologiqucs. Thesis, November 2006. 

[18] C. Poignard, P. Dular, L. Krahenbiihl, L. Nicolas, and M Schatzman. Methodes asymp- 
totiques pour le calcul de champs. In 5eme Conference Europeenne sur les Methodes 
Numeriques en Electromagnetisme (NUMELEC 2006), 29-30 Nov. 1 Dec. 2006. Lille, 
FRANCE. 

Centre de Mathematiques Appliquees, UMR CNRS 7641 and Ecole Polytechnique, 
91128 Palaiseau, poignard@cmapx.polytechnique.fr. 



